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Abstract 



Leaves of vascular plants are arranged regularly around stems, a phe- 
nomenon known as phyllotaxis. A constant angle between two successive 
leaves is called divergence angle. On the one side, the divergence angle «o 
of an initial pattern of leaf primordia at a shoot apex is most commonly an 
irrational number of about 137.5 degrees, called limit divergence. On the 
other side, the divergence a of a final pattern of leaf traces in the vascu- 
lar system of a mature stem is expressed in terms of a sequence of rational 



numbers, |, |, |, |, ^, called phyllotactic fractions. The mathematical 
relationship between the initial divergence ao, the final divergence a, and the 
number of internodes traversed by the leaf traces n c is investigated by means 
of a theoretical model of vascular phyllotaxis. It is shown that continuous 
changes of the trace length n c induce transitions between the fractional or- 
ders in the vascular structure. The vascular phyllotaxis transition suggests 
an evolutionary mechanism for the phenomenon of phyllotaxis. To provide 
supporting evidence for the model and mechanism, available experimental re- 
sults for fossil remains of Lepidodendron and the vascular structure of Linum 
and Populus are analyzed with the model. 
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Figure 1: (a) A typical pattern of leaf primordia (points) on a shoot apex with the initial 
divergence of «o = 1/ (1 +t) , or 360ao — 137.5 in degrees. The irrational number r ~ 1.618 
is the golden ratio defined by the proportion equation 1 : r = r — 1 : 1. The primordia are 
numbered in the reverse order of production. A solid spiral connecting all the primordia 
in the numerical order is the genetic spiral. Three dashed spirals (clockwise inward) 
and five dotted spirals (counterclockwise inward) are 3 and 5 parastichics, respectively. 
This pattern has a parastichy pair (3,5). (b) A typical pattern of leaves on a mature 
stem characterized with a divergence fraction of a = § (360a = 144°). Oblique strands 
diverging to leaves 1 and 6 are leaf traces. A solid spiral surrounding the stem is the 
genetic spiral. 



1. Introduction 



1.1. Review, background, and motivation 

Astonishing regularity manifested in plant architecture has fascinated var- 
ious fields of scientists for centuries. The regular arrangement of leaves, 
flowers and floral organs of higher plants is called phyllotaxis. A constant 
angle of rotation between two successive organs is called divergence angle, on 
which two apparently irreconcilable concepts have been in general use since 
the incept i on of quanti tativ e investiga ti ons o n phyllotaxis. 



Braunl ( I183ll . Il835l ) and ISchimperl (118351 ) noticed that divergence angle 



is various but not arbitrary. It is a fraction, or a rational number, a num- 
ber that can be expressed as the quotient — of two integers n and m. The 
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most widespread is the helical phyllotaxis, also called spiral or alternate phyl- 
lotaxis, in which stems bear a leaf per node. In the helical phyllotaxis, the 
numerator n and denominator m of the fraction normally are two alternate 
terms of a Fibonacci sequence, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, ■ • • . It is gen- 
erated by the Fibonacci recurrence relation that each number after the first 
two terms is the sum of the previous two numbers. The phyllotactic fractions 
I' I' I' I' 13 ' 2i' M' ' " com P r i se what is called the main sequence of phyl- 



lotaxis. A jr phyllotaxis is schematically shown in Fig. 1(b) In multijugate, 
verticillate or whorled phyllotaxis, where more than two leaves are borne at 
each node, the divergence angle is divided by the number of leaves in a whorl. 
In general, a plant stem is partitioned into nodes and internodes. A node is 
a point at which a leaf or leaves are attached, and an internode is a section of 
the stem between two successive nodes. In the fractional phyllotaxis, there 
are leaves aligned vertically above each other along a stem, as represented by 



leaves 1 and 6 in Fig. 1(b) A straight line connecting the superposed leaves 
is called an orthostichy. In the helical phyllotaxis, the denominator of the 
phyllotactic fraction is equal to the number of orthostichies. It is also the 
number of internodes between two adjacent leaves on an orthostichy. Thus, 



the 4 phyllotaxis in Fig. 1(b) has five orthostichies, 1-6, 2-7, 3-8, 4-9 and 



5-10, and five internodes separate leaves on each orthostichy. An imaginary 
spiral connecting all the leaves in the order of production is called the ge- 
netic, fundamental, generative, or ontogenetic spiral. The numerator of the 
fraction refers to the number of turns of the genetic spiral between the two 



adjacent leaves on an orthostichy. In Fig. 1(b), a solid spiral is the genetic 



spiral. From the leaf 6 to 1, the genetic spiral winds around the stem twice, 
the number two being the numerator of |. As remarked below, the phyl- 
lotactic fraction does not lose its significance even though vertical alignment 
is actually not exact but approximat e. 



In contrast, iBravais and Bravaid (118371 ) suggested that divergence angle 



is uniquely and invariably given by an irrational number, that is, a number 
which cannot be expressed as a fraction. The most typical angle of 360/(l+r) 
degrees is called the golden angle, where the irrational number r, called the 
golden ratio, golden mean, golden section, or extreme and mean ratio, is 
defined by the proportional relation 1 : r = r — 1:1. As the positive 
solution of the quadratic equation r(r — 1) = 1, it is given by 

t = n + 1 ~ 1.61803399 • • • . (1) 
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The defining equation is transformed to r _1 = 1/(1 + r _1 ). Recursive sub- 
stitution of t~ 1 in the left-hand side to the right-hand side gives an infinite 
continued fraction representation, 



1 



1 + t 



1 + 



The golden angle 360/(1 + r) = 360/r 2 is approximately 137.50776 degrees. 
By definition, the golden angle is the smaller angle created by sectioning the 
circumference of a circle (360 degrees) according to the golden ratio 1 : r, the 
golden section. A phyllotactic pattern with divergence equal to the golden 



angle is shown in Fig. 1(a) The ratio of the angle subtended by 1 and 2 to 
the angle between 1 and 3 is r, or Z102 : Z103 = 1 : r — 1 = r : 1, where O is 
the origin. Similarly, Z104 : Z204 = Z109 : Z409 = 1 : r, and so on. Thus, 
phyllotactic patterns with constant divergence equal to the golden angle have 
harmonious proportions. Patterns with an irrational divergence angle have 
no orthostichy in a strict sense, as no two leaves align vertically or radially. 
Instead, therefore, attention is directed to secondary spirals connecting posi- 
tionally nearby leaves, called parastichies. Like an orthostichy, a parastichy 



is characterized by a difference in number of leaves on it. In Fig. 1(a), the 
genetic spiral, three parastichies and five parastichies are drawn with a solid 
curve, dashed curves and dotted curves, respectively. Each of three paras- 
tichies 1-4-7-10, 2-5-8 and 3-6-9, is called a 3-parastichy. Hence there are 



three 3-parastichies and five 5-parastichies in Fig. 1(a), and the pattern in 



Fig. 1(a) is said to have a parastichy pair of (3,5), which is also denoted 
as (3 + 5) or 3 : 5. As a remarkable fact, parastichy numbers are almost 
always given by Fibonacci numbers. This is a mathematical consequence of 
the fact that divergence angle is almost always the special irrational number, 
the golden angle. The golden angle is also called the Fibonacci angle, for it 
is the limit angle of divergence for the phyllotactic fractions belonging to the 
main sequence; 

360 x | = 144, 

360 x § = 135, 
360 x | ~ 138.46, 
360 x £ ~ 137.14. 

The rational angles beyond ^ are practically indistinguishable from the 
'ideal' irrational angle of 137.507764- • • degrees. Therefore, it is argued that 
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what appear to be different rational angles are nothing but a single irrational 
angle disturbed by inevitable random errors. 

The seemingly conflicting views on the divergence angle, whether rational 
numbers or an irrational number, are a source of inspiration and confusion. 
In effect, they are not only compatible but both indispensable. On the one 
hand, the irrational number applies to the divergence angle of phyllotactic 
patter n s of und i fferentiated tissues at shoot tips or apical meristems (I Church 
(119041 ); iHirmerl (119221 . Il93ll )). Let us call it the initial divergence angle. It 
is commonly referred to as the ideal or limit divergence angle for the reason 
mentioned above. On the other hand, the rational (fractional) divergence 
applie s to phyllotax i s of le a ves, or primary vascular architecture on a mature 
stem (ILestiboudoisI (118481 ); iNagelil (118581 )). In the literature, the majority 
of studies discuss the former, i.e., the process of organ initiation, position- 
ing of the leaf primordia from which leaves will develop, and transitions of 
patterns at the shoot apical meristem. In recent years, substantial progress 
has been made in understanding plant hormonal factors that influence or 
control th e formation of lea f primordia and thei r arrangement on the apical 
meristem (IReinhardtl (120051 ); iKuhlemeierl (120071 )). In striking contrast, the 
fractional phyllotaxis of the mature stem have received less scholarly atten- 
tion, unfortunately. This is not because the latter is less important than the 
former. As a matter of fact, experimental findings on the close relationship 
between phyllotactic fraction and vascular organization ha ve been accumu - 
lated without being theorized from a general perspe c tive (ISterlingj (Il945h ; 
Girolamil (Il953h:ljensenl (fl9 68) ;lNamboodiri and Beckl (ll968l ); lL arson! (Il977l ); 



Beck et al.l (119821 ): iKirchofll (119841)) . 



Since the influential text by iHofmeisterl (118681 ). research into causal or 
dynamical mechanisms of primordia initiation has been the central pillar in 
the study of phyllotaxis. The empirical observation that new leaf primordia 
arise in the largest space between the older primordia is called Hofmeister's 
rule. What was originally a rule of thumb of botan i sts ha s been refined 
and developed into causal or dynamical models. lAiryl (118731 ) speculated on 
a causal mechanism in terms of geometrical objects in mechanical action. 



Schwendenerl (118781 ) put a similar idea on a more solid mathematical basis 



by regarding leaves on a stem as solid disks contiguously covering a cylinder 
surface of infinite length (Fig. [2]). In Schwendener's model, contiguous circles 
of a constant radius are arranged in a periodic pattern characterized with a 
given set of contact parastichy numbers. Then it is a purely geometrical 
problem to derive various mathematical relations for the divergence angle, 
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Figure 2: Schwendener's causal model. Right: Contiguous circles with decreasing radius 
are stacked on an unrolled surface of a stem cylinder. Contact parastichy numbers, or 
differences in the numbers of the circles in contact, change from (1,1) at the bottom to (5,8) 
at the top. Left: A mathematical relation between the divergence angle (the horizontal 
axis) and the radius of the contiguous circles (the vertical axis) is indicated with a solid 
zigzag curve starting from the top left corner (divergence of 180°, corresponding to the 
bottom part of the right figure) down to the golden angle 137.5° (the top of the right 
figure). The zigzagging is due to shifts in the contact parastichy numbers from (1,1) 
through (1,2), (2,3), (3,5), (5,8), (8,13), (1 3,21) to (21, 3 4). T he top branch for (1,2) 
extends from 180° to 128°34'. Adapted from ISchwendenerl f|l883h . 
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the radius of the circles, the girth of the cylinder, and the set of parastichy 
numbers. The radius of the circles is regarded as an independent variable, 
or a control parameter of the model. By letting the radius change continu- 
ously along the stem cylinder, the divergence angle varies concomitantly with 
the contact parastichy numbers according to the mathematical relations. As 
a remarkable result, the divergence angle converges toward the golden angle 
137.5° by decreasing the radius sufficiently slowly from a large initial value to 
a small constant value. Before attaining to the golden angle, the model pre- 
dicts that the divergence angle oscillates with decreasing amplitude (Fig. [2]). 
The decrease in the radius corresponds to decrease in relative size of leaf pri- 
mordia on the stem or apex. This is a brief summary of Schwendener's causal 
mechanism for the golden angle. The model is referred to as a mechanical or 
causal model of phyllotaxis. 



Related causa l models were discussed in depth by IDelpinol ( 118831 ) and 
van Itersonl (119071 ) . In recent decades, models of Schwendener and van Iter- 



son h a ve been elaborated on an d develop e d furt h er ma t hema ti cally (Adler 



(1974); iRothen and Kochl (I1989T): iLevitovl (119911): iKunzl (120011): lAtela et al 



(120021) ) or numerically dWilliams and Brittainl ( Il984h : Iflellwig et al.l (120061)1 



andeven realized dynamically in a physics laboratory experiment (IDouady and Couder 
( 119961 )). The causal models are founded on the basic assumption of causal 
determinism that a phyllotactic pattern is a result of causal interaction of 
pattern units. In particular, the position of an initiated leaf primordium is de- 
termined by the position of the older primordia according to (supposedly sim- 
ple) causal rules. The manner in which the units are arranged depends on the 
dynamic history of growth, or particularly on the course of changes in size of 
leaf primordia. Thus, a common key factor of the causal models is a gradual 
change in size of leaf primordia under mutual repulsion. Accordingly, there is 
a variety of causal models in which the repulsive interaction is ascribed not to 
the mechanical contact pressure as supp o sed b y Airy and Schwe n dener, but 
to a c h emical diffusion process (ISchoutd (I1913T) : iThqrnlevi (I1975T): iMitchison 



(1977) 



(1983) 



(1987) 



( 119941 ) 



Schwabe and Clewer ( 



"Veen and Lindenmayei (119771 ) : lYoung (1978) : iMarzec and Kappraff 



1984): I Chapman and Perrvl (119871): iRoberts 



Steeves and Sussex! (1989) ; lYotsumotol (119931 ): iKoch and Meinhardtl 



Meinhardt et al . (1998 )). There are anot h er causal models bas ed on 



physical (iHofmeisterl (Il868l): iGreen elaP (Il996l): iNewell et all (120081 )) and 



chemical ( Cummings and Strickland ( 1998| )) instabilities. Recently, more in- 
tricate models based on molecular-genetic experiments have been d iscussed 
dSmith et all (l2006al ]bh: Ijonsson et aD (120061 ): IShipman et~aD (l201l[ )). while 



geometrical mo d els have been used t o interpret patterns of real system s 
flMalvginl d2006h : iHotton et all (l2006h : IZagorska-Marek and Szpakl (l2008t )). 
All these causal models are based on the assumption that divergence angle 
is intrinsically variable and determined causally. 

In the recent literature, we have had few opportunities of finding phyl- 
lotactic fractions in use. Most theoretical and experimental works attach 
importance to parastichy numbers instead, and the phyllotactic fraction is 
not menti o ned o r regarded merely as an ap p roximation eve n if m e ntion ed 
flWilliamsl (1l974T ): ISteeves and Sussex] (119891 ): iLvndonl (Il990l ); Ijeanl (119941 ) I 
The fractional phyllotaxis is the original problem. There are some reasons 
for this trend. First, early researcher s did not apprecia t e the structural sig - 
nificanc e of th e phyllotactic fr a ction (IHofmeisterl (11868 ) : Ide Candolld ( 1l88ll ) ; 
Churchl (Il920h : Irlirmerl (Il922h : iRichardsl (Il95lh : ISnowl (Il955h l Second, the 



studies of vascular structure organization are comparatively so few in number 
that they are overshadowed by intensive research interests directed towards 
the shoot apical meristems. Third, Schwendener's causal model and its de- 
scendants are at variance with the fractional divergence. According to the 
model, the divergenc e angle varies depen ding on size of leaves, the vertical 
coordinate of Fig. [2j ISchwendenerl ( 118831 ) guessed that a fractional pattern 
would be made secondarily as a result of me chanic a l stra ightening of paras- 
tichous bundles connecting initiated leaves. iTeitzl (118881 ) confirmed indeed 
that the fractional phyllotaxis is accomplished by secondary torsion occur- 
ring in the vascular system during growth of the stem. When there is little 
or no secondary distortion for lack of subsequent growth or internodal elon- 
gation, the original pattern at the apex may grow to a similar pattern of 
mature organs. This holds true for the most eye-catching patterns of closely 
packed reproductive organs, which, therefore, are often compared favorably 
with numerically simulated outputs of causal models. Even then, the basic 
concept of the phy llotactic fraction may rema in significant internally in vas- 
cular connections ( Watson and Casper (jl984l )). In fact, a stem with short 
internodes takes a high-order fraction, which can be indistinguishable from 
the limit divergence of the undistorted stem. Thus, the phyllotactic fraction 
may not be judged by the external appearance alone. 

In vascular plants, each leaf is connected to the main stem vascular sys- 
tem through a strand of fluid-carrying vascular tissue called a leaf trace. 
A leaf may have several to many leaf traces. Leaf traces diverge from the 
stem vascular system some distance below or very near t he nodes at which 
they enter the leaves dLestiboudoisI (Il848f ): iNaeelil (Il858f ): iBeckl (l2010h ). At 
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the level of the shoot apex, leaf traces form parastichous strands winding 
obliquely round the stem axis. As the stem elongates, the leaf traces align 
up along the stem to make orthostichous bundles by forcing the whole stem 
to twist slightly from the original pattern (Fig. 1(b) ), thereby a phyllotactic 
pattern characterized by a phyllotactic fraction is established. In the final 
pattern, there is a definite relationship between the denomi nator of the phy l- 
lotactic fraction and the number of vascular orthostichies (Kirchofj f 1984T )). 
The term orthostichy might be misleading, because the straightened bun- 
dles still may maintain their tilted course. The fraction neatly represents 
geometrical arrangement of leaf traces, and the fractional order need not 
mean that leaves are positioned exactly vertically. Developmental sequences 
in differentiation and vascularization of leaf prim ordia are numerically corre - 
lated with the phyllotactic fr action of the shoot (jPriestley and Scottl (119361 ); 
Girolamil (119531 ); lEsaul ( 119651 )). Seemingly irregular rhythmical variations in 



various lengths of the external structure of a mature plant may be under- 
stood as a con se quence of a hidden phyllota c tic ord er in the vascular system 
( Unruh ( 1950 ); Kumazawa and Kumazawa ( 197ll )). There is evidence for 
restricted p athways of translocat i on of photosynthetic assimilates related to 
phyllotaxis ( Watson and Casper ( 19841 )). The patterns of translocation are 
sectorial, or the phyllotactic fraction has biological significance. The ob- 
servation most pertinent to the present work is the significant correlation 
existing between the phyllotactic fraction and the number of internodes tra- 
versed by leaf traces: The highe r phyl lo tactic fracti ons are associated with 
the longer leaf traces ( Girolami ( 1953 ); Esau ( 19651 )). While leaf traces of 
plants with helical phyllotaxis typically traverse more than one internode, in 
distichous phyllotaxis, a | phyllotaxis of two-ranked leaf arrangement, and 
i n verticillate phyl lotaxis, leaf traces are approximately one internode or less 
( IBeck et al.l ( 119821 )). Accordingly, low-order systems of a | and | phyllotaxis 
are seen on the stems of plants with long internodes, while plants with short 
internodes show high-order fractions, as remarked above. 

The stem vascular bundles, or axial bundles, and associated leaf t races 
compr ise sympodia, on the nature of which there are two perspectives ((Beck 
( 120101 )). In one vie w, the sympodia a re of cauline origin, or derived from 



(Il982h1. 



stem vascular tissue (IBeck et al 
origin, or derived from leaf traces (lEsaul (Il965l )) 



I n the o ther view, they are of foliar 
There are two different 



views on the causal relation between initiation of primordia and develop- 
ment of leaf traces or procambial strands, the strands differentiating into 
vascular bundles of xylem and phloem. In one view, the initiation of leaf 
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primordia brings about the differentiation of the leaf traces. Hence the ini- 
tiation of the leaf traces occurs basipetally, or in the direction from the leaf 
primordia toward the vascular system of the stem. In the opposing view, 
phyllotaxis of leaf primordia is dictated by the vascular org anization that 
has b een established before the leaf primordia are initiated (jLarsonl (119771 . 
19831 ) ). The latter is consistent with the observation that an incipient leaf 



trace d evelops acro p etally, or in the dire c tion t o ward the leaf primordium it 
serves ( Esaul ( 1965 ); Nelson and Dengler ( 1997t )). Priestley and Scott ( 1933 ) 
was criticized by ISnow and Snowl ( 11934 ) . They both do not cast doubt on 
Hofmeister's rule, that is, they share the causal view that phyllotaxis is a nat- 
ural consequence of growth and development of an individual plant. They 
differ in what they regard as a basic unit of phyllotaxis. The former adopts 
growth units including leaf traces, while the latter places primary emphasis 
on leaf primordia at the apex. Accordingly, the former and the latter at- 
tach little importance to the irrational and rational divergence, respectively. 
Thus, the causal view has been the paradigm of phyllotaxis. 

On the whole, causal models are successful in deriving indefinitely contin- 
uing stable systems, resembling actual phyllotactic patterns. From a compu- 
tational point of view, they are particularly appealing in that they provide 
us with programmable protocols leading to the golden angle. Irrespective 
of detailed mechanisms, however, realistic phyllotactic patterns are derive d 
based on the following observational fac ts (jVogell (119791 ) ; iRivier et al.l ( 119841 ) ; 
Prusinkiewicz and Lindenmayerl ( 1199 ll )): (i) Divergence angle is constant, 
(ii) The constant is the golden angle. For the sake of argument, the former 
is often taken so broadly that the constant may take any value. On this 
premise, phyllotaxis is rendered to a geometrical playground of mathemat- 
ics. There a re mathemat i cal ar g uments for (ii) based on the generali z ed hy - 
pothesis (D ride Candolld dl88lf ): ICoxet"en (Il972h : [E^ighl (jl972h : lRidleyl ( ll982f ); 
Marzec and Kapprafll <jl983h ). It is often stated in this regard that the golden 



angle is a special angle at which optimal packing is achieved. As a matter of 
fact, this is not true literally, for it is only under the constraint (i) that the 
golden angle may be said optimal and there is no a priori reason for the con- 
stancy. For living organisms, the property (i) is far from obvious and no less 
astounding than (ii), especially because the angular regularity may persist 
in spite of temporal irregularity. A time interval between the formation of 
successive leaves is called a plastochron, which is used as a morphological or 
developmental time scale. Plants grown in different environmental conditions 
may be compared in plastochron units but not in physical time units. The 
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fact that unit of time is a plastochron and duration of a plastochron is not 
constant in physical time poses a problem for realistic causal models based 
on physical time. 

Phyllotactic patterns at a shoot apex are more regular than those on a 
mature stem, because internodes tend to be elongated less regularly on the 
mature stem. As a matter of fact, the exact level in the stem at which a bifur- 
cation or recombination of vascul a r bun dles takes place is not an important 
morphological constant ( iDormerl ( 119721 ) ). Accordingly, trace lengths may 
vary arbitrarily along the stem. For this reason, it is often argued t hat one 
should devote onese lf exclusively to the study of the growing apex (j Church 
(1l904T ): ISnowl ril955j )). Nonetheless, further mathematical relations for the 
spiral patterns at the apex can be derived by assuming a stronger mathemat- 
ical constraint of exponential growth, according to which the leaf primordia 
are arra n ged o n logarithmic spi r als in a centric r epres e ntatio n (Fig. 1 (a) ) 



flChurchl (119041 ): iRichardsl fll95l|); iThomad (1 19 75( 1: Ueanl (|1994j)). In the ex- 
ponential growth, the ratio of the distances from the center of the apex to 
two successively numbered primordia is a c onstant, called the plastochron 
ratio. For a fixed value of divergence angle, IRichardsl (ll95l[ ) has advocated 
the use of a phyllotaxis index defined in terms of the plastochron ratio (cf. 
( lA.ip ). The index is used to designate two sets of parastichies intersect- 
ing orthogonally. For instance, for Fig. 1(a), the plastochron ratio is 1.2, 
the phyllotaxis index is 3, and (3, 5) parastichies cross at right angles. A 
fractional value of the index, such as 2.7, means that no two parastichies 
are orthogonal. In this geometrical model, a shift in parastichy numbers, 
e.g. from (3,5) to (5,8), a phenomena called rising phyllotaxis, is related 
to a variation of the plastochron ratio, or the exponential growth rate. The 
model has been g eneralized to allow for other cons t ant divergen ce angles than 
the golden angle ( Richards ( 1951 ): Thomas ( 1975 ): Jean ( 19941 )). In contrast 
to these geometrical models based on constant divergence angle, there ex- 
ist geometrical causal models in line with Schwendener's model, which aim 
at deriving the limit divergence angle by assuming variable divergence an- 
gles depending on plastochron, the plastochron rat i o, and their own rules 
(Ivan Itersonl (119071 ) : IWilliamd (119741 ) : lEricksonl (119831 ) : IWilliams and Brittain 



(119841 ) ) . For a vegetative shoot, a plastochron index is defined in terms of 



lengt h of leaves, and a leaf on a sh oot is labeled with a leaf plastochron in- 
dex ( Erickson and Michelini ( 19571 )). A developmental index of this kind is 
i ndisp ensable for the systematic study of plant development (IMeicenheimer 
(120061 )). The exponential growth is a practically useful approximation in 
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dealing with young organs and early stages of development, although it is 
neither essential nor peculiar to phyllotaxis. 

Despite the apparent success of causal models, neither their intrinsic 
mechanisms nor predictions have yet been subjected to experimental tests 
specifically. To name several problems on a descriptive level, existing causal 
models that explain all types of observed patterns cannot help predicting also 
a multiplicity of unreal or too rare patterns. Even when they are capable 
of deriving normal patterns, they are not free from instabilities apparently 
irrelevant to living organs. Causal models in general are confronted with 
a subtle trade-off. Normal phyllotactic patterns must be stable enough to 
account for the current prevalence in nature, while they cannot be quite sta- 
ble in order to allow for many other exceptional ideal angles just to such 
a degree that they are actually existent. In short, rare patterns should be 
neither too common nor too rare. It is not clear how and why this subtle 
balance between stability and instability is maintained universally, since fine 
control of relative size of phyllotactic units depends not only on species but 
individual plants or even on parts of the individual plant. We get puzzled 
all the more by the observations of more frequent occurrence of rare patterns 
among fossil plants. 

Causal models, whether physical or chemical, provide dynamical schemes 
of self-adjusting the system under the influence of the older leaf primordia. 
On the premise that divergence angles between successive leaves are freely 
variable by nature, they aim to derive a special angle, normally the golden 
angle, toward which the variable divergence angles tend ultimately. They do 
not assume any special constant divergence a priori. For the very reasons, 
they are likely to be beset with a fundamental difficulty in protecting the sys- 
tem against disturbance. In this regard, Hofmeister's empirical rule is often 
overestimated. Observed patterns satisfy Hofmeister's rule, but Hofmeister's 
rule is not sufficient for observed patterns. Hofmeister's ru le does not imply 
the periodic appearance of new primordia (iKirchofil ( 120031 )). nor does it en- 
sure precise regulation of the divergence of 137.5° (cf. Fig. 12(a)). It is not 
difficult to draw an unreal pattern according to Hofmeister's rule. The re- 
markable empirical fact is rather that divergence angle during steady growth 
seems always regulated stably to one of special angles closely related to the 
golden ratio. In fact, if a causal rule is to be strictly applied throughout, 
fluctuations in size of the domain of influence of a leaf primordium should in- 
evitably leave behind everlastin g; irregularities propagated in the developing 
pattern (ISnow and Snow! (119621 )). Mathematically, the instability is a general 
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consequence of the fact that the number of possible phyllotactic configura- 
tions proliferates as relative size of phyllotactic units decreases. According 
to causal interpretations, higher phyllotaxis becomes more vulnerable. The 
difficulty may not be obvious if one were interested only in low-order patterns 
like a (2,3) and (3, 5) phyllotaxis, but it should become conspicuous when 
dealing with a higher order pattern which requires higher precision mainte- 
nance. Besides this stability problem for high-order patterns, causal models 
have another difficulty for low-order patterns (Sec. |3]). 

There are apparent geometric correlati ons between p a rastichy number s 
and rela t ive siz e of primordia on the apex (I Church! ( 11904] ); iRichardsl (119511); 
Kirchofll (|2003j )) and between leaf arcs and the plastochron ratio (IRutishauserl 



( 119981 )). Causal models implement them as causal relationships with the in- 
tention of proving that a phyllotactic pattern, especially the divergence angle 
of 137.5°, is a necessary consequence of changes in the causal agent, relative 
size of leaf primordia. According to this interpretation, divergence angles 
and contact parastichies must depend not only on the shape of primordia 
but on the geometry of the su rface on wh i ch the y are located. The depen- 
dence has been investigated by Ivan Itersonl (119071 ) on the assumption that all 
the primordia keep a common shape while they are allowed to change their 
sizes. So far, however, no direct evidence has been provided to support the 
presumed causal relationship. As a matter of fact, there are very few studies 
in which sufficiently detailed data are obtained to ma ke a close compari- 
son with the models possible or useful ( lEricksonl (119831 ) ) . In particular, the 
prediction of causal models that rising phyllotaxis, or change in parastichy 
numbers, should accompany wide variations and abrupt turns of divergence 
angle, as indicated in Fig. |2l has not been supported unequivocally. On the 
contrary, the success of Richards' model indicates the exponential growth 
with consta nt diverg e nce a ngle irrespective of whether parastichy numbers 



rise or fall. I Church! (119041 ) refuted Schwendener's model by counterexam- 



ples showing normal spiral patterns of circular primordia whose positions are 
widely separated. In comparing treated plants, iMaksymowych and Erickson 
(119771 ) found no significant change in divergence angle in a correlation dia- 
gram for_the plastochron ratio and divergence angle. Statistical analysis of 



Fujital (119391 ) has revealed that divergence angles do not depend so much on 
parastichy numbers as expected from causal models (Jean] (119861 )). There is 
clear evidence against the basic assumption that the primordia size is the 
causal factor of divergence angle. A plant appears to accomplish geometri- 
cal correlations in a phyllotactic pattern by adapting the size and shape of 
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l eafy o rgans as if it knows the end pattern at which it aims. ISnow and Snow 



(119621 ) observe that the secondary extension of a leaf base adjusts itself so 
that divergence angle is little affected in spite of artificial disturbances. This 
observation, despite the authors' claim, undermines their space-filling mech- 
anism that the leaf base extension regulates the phyllotactic pattern. To the 
contrary, apparent causal changes in the position, size and shape of leaves 



( Church 


f 


1904); 


Richards 


( 


1948 





there are cases in which physical or chemical contact pressure may induce 
secondary displacement of compactly packed lateral organs. 

Natural selection plays no role in causal interpretations of phyllotaxis. If 
one supposes to the contrary that natural selection holds the key to under- 
standing the golden angle at the shoot apex, then one should investigate a 
special effect of the special angle, instead of its cause. In other words, one 
should look for the distal or ultimate cause of the special angle, instead of 
the proximity cause. This sort of theory intends to explain special traits 
not in terms of immediate physiological factors, but in terms of evolutionary 
forces acting on them. It aims at a full understanding of the phenomena 
at a phenomenological level, independently of whatever physiological mech- 
anisms may be involved. There is a long history of investigations into selec- 
tive advantage of the observed divergences based on the external structure. 
It goes as follows: common phyllotactic patterns distribute leaves as evenly 
as possible and maximiz e expo sure of leaves to enhance the capacity to in- 
tercept sunlight ( Wright ( 1873h ) . Such an argument is unpromising because 



leaves are aligned in vertical ranks. Indeed, changes in leaf sh ape an d stem 



lengt h can compensate for the negative effects of leaf overlap (jNiklasl ( 119881 . 



19981 )). For this obvious reason, it is often argued to the contrary in favor 
of the 'most irrational' divergence angle; no two leaves lie pre cisely under 
one another when divergenc e angle is equal t o the golden a ngle (de Candolld 



(ll88lh : IWies^eil (Il875l. Il907h : ICoxeteri fll972h : iLeighl (Il972|) : iTakenakal (I l994 h 



Pearcv and Yand (ll998f ): IValladares and Britesl (120041 ); King et all (120041 ); 
Brvntsevl (120041)'). There is also a long history of criticism of this view 
( Thompson ( 1917h ) . In the first place, the golden angle is not a general 
rule for mature shoots, and the light-capture mechanism deepens the riddle 
of the common occurrence of a | phyllotaxis. In general, existing theories 
tend to argue for the uses of irrational angles without regard to the uses of 
rational angles or vice versa. 
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1.2. Aim and scope of this paper 

Let us direct attention to the internal structure, the vascular system. 
Mathematical interrelationship between the initial (apical) and the mature 
(vascular) phyllotactic pattern seems to have not been discussed experimen- 
tally nor theoretically. This paper develops a theory of vascular phyllotaxis 
to fill in the gap between the two distinct but intimately related phenom- 
ena. A physical mo del has been described mathematically in the previous 
paper dOkabd fl201lh . However, the model was abstract and its relevance 
to real phenomena was not clearly elucidated. The aim of this paper is to 
develop the model to show its experimental validity and relevance. This is 
relevant to a fundamental problem of phyllotaxis: Is phyllotaxis determined 
causally or genetically? In contrast to numerous models holding the causal 
view, the present model is based on the genetic perspective that special num- 
bers in phyllotaxis are primarily of genetic origin, so that it is assumed that 
constant primordial divergence angle during steady growth is genetically de- 
termined. According to the model, the effect of constant divergence angle 
is investigated, and what value of the constant is advantageous is settled. 
This work is not concerned about transient fluctuations of divergence an- 
gle during ontogeny. Therefore, the model is compatible with any physical 
or chemical causal models for the positioning of leaf initiation at the shoot 
apex, although the limit divergence angle at the apex is interpreted totally 
differently. The special angle is not an inevitable consequence of ontogenetic 
dynamics, whether physical or chemical. It is regarded as a heritable trait of 
a plant. It is supposed that once there was a wide variation in the traits of 
individuals, or there have formerly been wide variations of divergence angles. 
The special limit divergences found in nature have survived natural selection. 
This conforms with the traditional view that biological features that are un- 
der tight genetic control and that hay e very narrow ranges of variation are 
believed to be adaptive (INiklas dliiazl)). Although the author believes that 
the premise of the model, divergence angle as a trait of a plant, is not only 
plausible but supported by circumstantial evidence, it has not been unani- 
mously accepted at present. It may be verified or refuted experimentally in 
the future. 

For the efficient transport of materials throughout an indefinite number 
of leaves attached to a stem of a finite cross section, the leaves should be 
aligned along a finite number of 'orthostichious' bundles. At this point, a 
whole number enters the theory. There are modes of orthostichous order de- 
pending on the initial arrangement and length of leaf traces. The number of 
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vascular orthostichies may be increased or decreased, but not arbitrarily. By 
regarding a leaf primordium and the leaf trace (s) associated with it as the 
fundamental unit of phyllotactic patterns, a mathematical correspondence 
is derived between the divergence angle of the initial phyllotactic pattern, 
360«o degrees, and the phyllotactic fraction a of a mature pattern, where 
the number of internodes traversed by the leaf traces, n c , plays a pivotal 
role. As a general rule, it has been known that ph yllotact i c frac ti on of a 
vascul a r plant may v a ry sequent i ally d u ring growth (|Braunl (118351 ) ; ISkutch 
fll927h : lAllardl (Il942h : IPulawskal (Il965h : iLarsonl (Il977h l By means of the 
mathematical relation between «o (an irrational number) and a (rational 
numbers), it is shown that changes in n c cause the phyllotactic transitions 
in a. As a natural consequence, an evolutionary mechanism for the phe- 
nomenon of phyllotaxis is suggested. Supporting evidence for the model and 
the evolutionary mechanism is presented by analyzing experimental results. 

In Sec. [2l a model and results used in the following sections are presented 
by means of figures and tables without using mathematics. Tables [TMT81 have 
not been presented before. 

In Sec. 0, observed precision of the initial divergence a is explained by 
means of a correlation predicted between the range of «o and the highest- 
order fraction a. In short, divergence angle «o of a system with a high 
phyllotactic fraction a should be accurately controlled in order to avoid un- 
necessary changes in vascular structure. 

In Sec. IU phyllotaxis of Lepidodendron by iDicksonl (118711 ) is analyzed. 
Diversity of phyllotaxis is discussed as a result of ineffective selective pres- 
sures. 

In Sec. 0, the vascular structure of Linum usitatissimum by iGirolami 
(119531 ) is investigated. Various relations between phyllotactic fraction and 
parastichy numbers, the phyllotactic fraction a and the length per internode 
of leaf traces n c , and directions of parastichies and the genetic spiral are 
pointed out. 

In Sec. El the phyllotactic transition of Populus deltoides by ILarsonl (119771 ) 
is analyzed. It is shown that a continuous change in length of leaf traces 
causes the discontinuous effect of the phyllotactic transition in the vascular 
structure. 

In the appendix, a relation between the trace length n c and the plas- 
tochron ratio a is discussed to indicate that the former serves as a useful 
developmental index for the mature stem as the latter is used for the apex. 
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Figure 3: Phyllotactic patterns of leaf traces with a length of n c = 4 before and after 
secondary torsion are arranged side by side. A dotted and dashed line of each figure 
represent a vertical cut of a cylinder surface unrolled, (a) A pattern with initial divergence 
of 36Ctao — 137.5° (eta = l/(l+r) ~ 0.382). (b) The final pattern of a fractional divergence 
a = | resulting from (a). Leaf traces in the upper part move rightward while the pattern 
(a) becomes (b), thereby five 5-parastichies in (a), such as 1-6-11-16-21, align themselves 
to make five orthostichies in (b). 



2. Model 

A regular helical pattern of leaf traces is schematically plotted as a lattice 
of line segments on an unrolled surface of a cylinder. The divergence angle 
of the initial pattern is denoted as 360«o in degrees, which is assumed to be 
less than 180 degrees, i.e., < a < | without loss of generality. Fig. 3(a) 
presents a typical pattern for 360«o — 137.5° («o — 0.382). The length of 
leaf traces measured in internodes is denoted as n r in accordance with the 



previous notation ( lOkabd ( ]201ll )). As in Fig. 1(b) n c = 4 in Fig. [3J The 
trace length n c need not be an integer; n c is the average number of leaf traces 
cut by a transverse section (Fig. HJ). As the number in a section is an integer, 
this method gives a good estimate of n c particularly for n c ^> 1. The model 
comprises two parameters «o and n c . For the sake of argument, patterns with 
constant values of them are considered below. Effects of their fluctuations 
may be discussed based on results to be obtained. 

The leaf traces repel with each other laterally to arrange themselves in 
an orthostichous pattern. Th e mut u al int eraction is likely to be regulated 



by the plant hormone auxin (jBecki ( 120101 )). Fig. 3(b) is the final pattern 
resulting from Fig. 3(a) Divergence of the final pattern is expressed in terms 



of the phyllotactic fraction a. The pattern of Fig. 3(b) is characterized with 
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Figure 4: For three transverse sections of a pattern of leaf traces with a length of n c = 4.3, 
the number of the traces in each section is indicated on the right-hand side below the cut 
line. The number averaged over sections should approach n c . 



a 



g. In Fig. 3(a), there are five parastichies of 1-6-11-16-21, 2-7-12-17- 
22, 3-8-13-18-23, 4-9-14-19-24 and 5-10-15-20-25, each of which is called a 
5-parastichy. The five 5-parastichies are lined up vertically to make five 
orthostichies of the | phyllotaxis in Fig. 3(b) In the patterns of Fig. [31 the 



next visible parastichies are 3-parastichies (1-4-7-10-13-16-19-22-25, 2-5-8-11- 
14-17-20-23 and 3-6-9-12-15-18-21-24) and 2-parastichies (1-3-5-7-9-11-13-15- 
17-19-21-23-25 and 2-4-6-8-10-12-14-16-18-20-22-24). As these parastichies 
remain conspicuous in the two patterns, both patterns may be referred to 
as having a parastichy pair of (2,3). Thus, for n c = 4, there is one-to-one 
correspondence between ao ~ 0.382 (angle of 360ao — 137.5°) of the initial 
pattern and a = | of the final pattern. In a similar manner, a is obtained 
for arbitrary values of ap and n c . In deed, we get a — f insofar as 3 < n c < 5 
and \ < a < 



sec 



Okabd (120111 ) for the mathematical implementation). 
Below we discuss phyllotactic changes in a that occur when n c and ao are 
set out of their respective ranges. 

For a fixed value of «o — 0.382, Fig. 5(a) is for n c = 7 in comparison 
with Fig. 3(a) for n c = 4. As the traces of length longer than five internodes 
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Figure 5: Change in a phyllotactic pattern of leaf traces with a length of n c = 7 (cf. 
Fig. El), (a) The initial pattern with a = 1/(1 +r) = 1/(2 + (360a ~ 137.5°). (b) 
The final pattern with a = |. 



for Fig. 5(a) , while 



cannot be aligned in five orthostichies, we obtain a = 

a = | in Fig. 3(a) Thus, it is explained that a higher phyllotactic fraction 
is obtained for a longer length of leaf traces. Phyllotactic transition from 
a = | to a = | occurs when n c increases past a threshold value of n c = 5. 
Experimental evidence of this transition is presented below in Fig. [T7J 

For a fixed value of n c , the phyllotactic fraction a depends on the initial 
divergence ocq. For n c = 4, Fig. [6] is for «o — 0.276 (angle of 99.5°), which 
is compared with Fig. [3] for «o — 0.382 (137.5°). The former leads to a final 
pattern of a = | in Fig. 6(b), while the latter gives a = | in Fig. 3(b) 
In fact, there are three fractional patterns conceivable for n c = 4 
a = | for < a < |, (b) a = | for ~ < a < | an d (c) a 



a 



namely 
= | for 

5 



| < a < f . 

Every phyllotactic fraction for a has its own ranges of values for ao and n c . 
The mathematical correspondence is presented succinctly as a tree diagram 
in Fig. [3 For instance, Fig. [7] gives the conditions \ < a < | and 4 < n c < 7 
for a = |. For the former inequalities, the boundary fractions | and | lie 
below | in Fig. [71 The latter condition 4 < n c < 7 is reasoned from the 
vertical coordinate n c = 4, 5 and 6 of three |'s in Fig. [7J The phyllotactic 
sequence of fractions derived from an arbitrary value of initial divergence «o 
may be traced by climbing up the tree of Fig. [7] along the vertical line at «o- 
For a ~ 0.382 (137.5°), the main sequence ^, ~, |, |, ^, • • • is obtained 
in the increasing order of n c . The tree diagram extended for all values of n c 
includes all conceivable phyllotactic fractions. 
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Figure 6: Change in a phyllotactic pattern of leaf traces with a length of n c = 4 (cf. 
Fig. El), (a) The initial pattern with a Q = 1/(3 + r" 1 ) ~ 0.276 (360a ~ 99.5°). (b) The 
final pattern with a = j. 



For the sake of convenience, let us introduce shorthand notation for the 
irrational numbers found in nature, 



n 



1 



n + t 



n + 



1 



[n, m, I] 



m + t 1 
1 



(2) 



n + 



m + 



l + T 



-1 



and so on, where n, m, I are positive integers. With this notation, «o = [2] = 
1 / (2+r~ 1 ) = 1/ (1+t) gives the main sequence. The last equality holds by the 
definition of r in ([T|). Note that the pattern with a® = [1] (360ao = 222.5°) 
is nothing but the mirror image of «o = [2] (137.5°), because [1] = 1 — [2] 
or 222.5° = 360° — 137.5°. For future reference, Tables [TKfT8l are provided 
for the initial divergence ao given by typical irrational numbers. These are 
not exhaustive, but they include almost all phyllotactic fractions observed in 
nature. The main sequence is presented in Table [U In the second column for 
a = 1, '(2, 3)' in the second row represents the parastichy pair corresponding 



to the fraction 



5 ' 



in the third row abbreviates 3 < n r < 5, where 5 for 



the upper limit is taken from the next column, and \ ~ |' in the fourth row 
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Figure 7: Tree diagram for the phyllotactic fraction a. The horizontal axis is the initial 
divergence ao, and the vertical axis is the number of internodes traversed by leaf traces 
n c . Numbers in parentheses below each fraction are the parastichy pair corresponding 
to the fraction. By way of explanation, let us take I in the right-bottom quoter as an 
example. The fraction | with the parastichy pair (3,5) is in three different positions at 
n c = 5,6 and 7. By means of lower order fractions lying below them, the fraction | is 
bracketed between | and |. Therefore, we obtain a — | with the parastichy pair (3,5) 
insofar as i < ao < § and 5 < n c < 8 (Table [T]). Similarly, we hnd a — with (5, 7) for 
| < ao < I and 7 < n c < 12 (Table E]). 
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Table 1: Parastichy numbers (in parentheses) and ranges of n c and ao for the phyllotactic 
fractions a belonging to the main sequence with the limit divergence of ao = [2] = 1/(2 + 
r -1 ) ~ 0.3820 (360ao — 137.5°). The parastichy pairs are generated from the seed pair 
(1,2) for a = | (not shown) by a Fibonacci recurrence relation. The golden angle ao = [2] 
and the main sequence are called the Fibonacci angle and the Fibonacci sequence. 
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Table 2: Table for the limit divergence of a = [3] = l/(3+r- 1 ) ~ 0.2764 (360a - 99.5°). 
The parastichy pairs are generated from the seed pair (1,3) by a Fibonacci recurrence 
relation. The sequence 1,3,4,7,11,- • • is called the Lucas sequence or the first accessory 
sequence. 



indicates | < a < \. As an example, let us take a fraction a = ||. It is found 
in the eighth column of Table |2j from which the conditions 47 < n c < 76 and 
< a < || are read. These results are used in the next section (Table [2"0]) . 
The tables show that the denominator of a fraction a and the parastichy 
numbers are correlated with the threshold numbers for n c . The numbers 
comprise a characteristic sequence of integers. The main sequence in Table [1] 
is characterized with the Fibonacci sequence of 1, 2, 3, 5, 8, 13, 21, • • • , while 
Table [6] has a sequence of 1, 2, 3, 2, 5, 7, 12, 19, • • • , which is sometimes 
called the lateral sequence. Three consecutive numbers of a sequence satisfy 
the Fibonacci recurrence relation (2 + 5 = 7, 5 + 7 = 12, 7+ 12 = 19), except 
for the first several numbers (like 1, 2, 3 in the lateral sequence). Therefore, 
each phyllotactic sequence is referred to by a pair of seed integers for the 
Fibonacci recurrence formula. The seed pair of each table, such as (2, 5) in 
Table El is highlighted in boldface. 
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"7v 1 1 ? A A A 15 21 — 

" 4 5 9 14 23 37 60 97 

(1,3) (1,4) (4,5) (5,9) (9,14) (14,23) (23,37) (37,60) 
n c 3~ 4~ 5~ 9~ 14~ 23~ 37~ 60~ 

3 U _ 4 5 ~ 4 5_Q 14 _ 9 14 _ 23 37 _ 23 37 _ 60 



Table 3: Table for the limit divergence of a = [4] = l/(4+r- 1 ) ~ 0.2165 (360a ~ 78.0°). 
The parastichy pairs except (1,3) are generated from the seed pair (1,4) by a Fibonacci 
recurrence relation. The sequence 1,4,5,9,14,- • • is called the second accessory sequence. 



I I I A A A A A 

4 5 6 11 17 28 45 73 

(1,3) (1,4) (1,5) (5,6) (6,11) (11,17) (17,28) (28,45) 





3~ 4~ 5~ 


6~ 


11~ 17~ 


28~ 


45~ 


«o 


o~i ~ ^ ~ -| 


i i 

— rs_/ — 

fi 5 


12 3 2 
fi ~ 11 17 ~ 11 


3 5 
17 ~ 28 


8 5 
45 ~ 28 




Table 4: a = [5] = 


: 1/(5 + T 


- 1 ) ~ 0.1780 (360a - 


64.1°). 




a 


111 
(1^3) (1*4) (1?5) 


1 

(M) 


2 3 

13 20 

(6,7) (7,13) 


5 

33 

(13,20) 


8 
53 

(20,33) 




n c 


3~ 4~ 5~ 


6~ 


7~ 13~ 


20~ 


33- 


a 


o~i o ~ ^ 0~-^ 


~ - 

fi 


1112 

— r*^ — — r*^ 

7 fi 7 13 


A ^ A 

20 13 


A ^ A 

20 33 




Table 5: a = [6] = 


1/(6 + T~ 


- 1 ) ~ 0.15112 (360a ^ 54.4°). 




a 


2 3 A 
(2*3) (2^5) (5,7) 


8 
19 

(7,12) 


13 21 
31 50 

(12,19) (19,31) 


34 
81 

(31,50) 


55 
131 

(50,81) 




n c 


3~ 5~ 7~ 


12~ 


19~ 31~ 


SO- 


81~ 


a 


112 12 3 

— rs_/ — — rs_/ — — rs_/ — 


A ^ 3 


A ^ A 13 _ A 


ls ~ 21 


34 21 

r-u 


.3 2 5 2 5 7 


12 7 


12 19 ,31 19 


31 50 


81 50 



Table 6: a = [2,2] = 1/(2 + 1/(2 + t" 1 )) ~ 0.4198 (360a ^ 151.1°). The parastichy 
pairs except (2,3) are generated from the seed pair (2,5). The sequence 2,5,7,12,19,- • • is 
called the first lateral sequence. 



23 



(X 


l 


2 

(3^4) 


3 5 
10 17 

(3,7) (7,10) 


8 13 

27 44 

(10,17) (17,27) 


21 

(27,44) 


34 
115 
(44,71) 


a 


~ \ 


4~ 
i i 

4 3 


10~ 

2 12 3 
7 3 7 10 


17~ 27~ 

5 3 5 8 
17 10 17 27 


44~ 

13 8 
44 27 


71~ 

13 21 
44 71 




Table 7: a = 


[3, 2] = 1/(3 + 1/(2 + 


t- 1 )) ~ 0.2957 (360a ~ 106.4°) 




a 


2 


3 

(2^5) 


4 7 

9 16 

(2,7) (7,9) 


11 18 
25 41 

(9,16) (16,25) 


29 
66 

(25,41) 


47 

107 

(41,66) 


n c 
a 


3~ 
i i 

— r*~/ — 

3 2 


5~ 

2 1 

— r*^ — 

5 2 


3 13 4 

— rs_/ — — rs_< — 

7 2 7 9 


16~ 25~ 

X „ 4 X ^ 11 

16 9 16 25 


18 11 

rs_/ 

41 25 


66~ 

18 29 

r^u 

41 66 


Table 


8: ao = 


[2,3] = 


1/(2 + 1/(3 + T- 1 )) 


~ 0.4393 (360a ~ 


158.1°). The sequence 



2,7,9,16,- • • is called the second lateral sequence. 



~~ 2 IT JL X II IS 2TJ JTT 

" 7 10 13 23 36 59 95 154 

(3,4) (3,7) (3,10) (10,13) (13,23) (23,36) (36,59) (59,95) 



n c 
a 


1121 

4 3 7 3 


10~ 13~ 23- 

3 13 4 7 

r^-> — rs_/ rs_/ 

10 3 10 13 23 


36~ 59~ 95~ 

4 7 ^ 11 18 ^ 11 18 ^ 29 
13 23 ~ 36 59 ~ 36 59 ~ 95 




Table 9: a = 


[3, 3] = 1/(3 + 1/(3 + T" 1 ))* 


0.3052 (360a ~ 109.9°). 


a 


1 1 

(M) (M) 


1 A. A. 

9 13 22 

(4,5) (4,9) (9,13) 


_8_ 13 21 

35 57 92 

(13,22) (22,35) (35,57) 


n c 
a 


3~ 4~ 
~ I ~ - 

u 3 u 4 


5~ 9~ 13~ 
11212 3 

— r^j — — rsj — — rs_/ 

5 4 9 4 9 13 


22~ 35~ 57~ 

5^3 5^8 13 ^8 
22 13 22 35 57 35 




Table 10: a = 


= [4, 2] = 1/(4 + 1/(2 + r" 1 ))^ 


i 0.2282 (360a - 82.2°). 


a 


(2*3) (3*5) 


4 7 11 

11 19 30 

(3,8) (8,11) (11,19) 


18 29 47 

49 79 128 

(19,30) (30,49) (49,79) 


n c 
ao 


3~ 5~ 
1112 

— rsj — — rvj — 

3 2 3 5 


8~ 11~ 19~ 

1 3 4 3 4 7 

— rsj — rs_/ — rs_/ 

3 8 11 8 11 19 


30~ 49~ 79~ 

11 7 11 18 29 18 
30 ~ 19 30 ~ 49 79 ~ 49 



Table 11: a a = [2,1,2] = 1/(2 + 1/(1 + 1/(2 + T" 1 ))) ~ 0.3672 (360a ~ 132.2°). 
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a 


•2 

&) 


3 

(2^5) 


5 

12 

(5,7) 


7 

17 

(5,12) 


12 

29 

(12,17) 


19 

46 

(17,29) 


31 

75 

(29,46) 


50 

121 

(46,75) 


n c 
a 


3~ 
1 i 

— rsj — 

3 2 


5~ 

2 1 

— rsj — 

5 2 


7~ 

5 7 


12~ 
2 _ A 

5 12 


17~ 

7^5 
17 ~ 12 


29~ 

X _ 12 
17 29 


46 ~ 

19 ^ 12 
46 ~ 29 


75- 

19 _ 31 
46 75 



Table 12: a Q = [2,2,2] = 1/(2 + 1/(2 + 1/(2 + r" 1 ))) ~ 0.4133 (360a ~ 148.8°). 



~~ I 2 IT jT X TT TS W 

a 4 7 11 15 26 41 67 108 

(1,3) (3,4) (4,7) (4,11) (11,15) (15,26) (26,41) (41,67) 

n c 3~ 4~ 7~ 11~ 15~ 26~ 41~ 67~ 

1 1 1 2 1 3 4 3 4 7 11 7 11 18 

— r^j — — r^j — — r^j rsj rsj rs_/ rsj 

4 3 A Z A LI 15 1J 15 2fi 4J 2fi 4J S7_ 



ao 



Tabic 13: a = [3, 1, 2] = 1/(3 + 1/(1 + 1/(2 + t" 1 ))) ~ 0.2686 (360a ^ 96.7°). 



n I 1 A. A _L I? W jK. 

" 4 7 10 17 24 41 65 106 

(1,3) (3,4) (3,7) (7,10) (7,17) (17,24) (24,41) (41,65) 
n c 3~ 4~ 7~ 10~ 17~ 24~ 41~ 65~ 

nil 12 12 52 57 57 12 19 12 

"=0 u 3 4 3 7 3 7 17 7 17 24 17 24 41 65 41 

Table 14: a a = [3, 2, 2] = 1/(3 + 1/(2 + 1/(2 + r- 1 ))) ~ 0.2924 (360a ^ 105.3°). 



a 


•2 

(is) 


3 

(3^5) 


4 
11 

(3,8) 


5 

14 

(3,11) 


9 

25 

(11,14) 


14 

39 

(14,25) 


23 

64 

(25,39) 


37 

103 

(39,64) 


n c 
a 


3~ 

1 i 

— rsj — 

2 3 


5~ 

1 2 
— rsj — 

3 5 


8~ 

1 3 

— rsj — 

3 8 


11~ 

1 4 

— r*u 

3 11 


14~ 

5 4 

r\_/ 

14 11 


25~ 

5 9 

r\j 

14 25 


39~ 
M „ A 

39 25 


64~ 

M ^ 23 
39 64 



Table 15: a Q = [2,1,3] = 1/(2 + 1/(1 + 1/(3 + T" 1 ))) ~ 0.3593 (360a ~ 129.3°). 
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a 


2 


3 

(2^5) 


4 


7 

16 

(7,9) 


10 

23 

(7,16) 


17 

39 

(16,23) 


27 

62 

(23,39) 


44 
101 
(39,62) 




(2*3) 


9 

(2,7) 


n c 
a 


3~ 
i i 

3 2 


5~ 

2 1 

5 2 


7~ 

3 1 
7 2 


9~ 

3 4 
7 9 


16~ 

3^7 
7 16 


23~ 

10 ^ 7 
23 16 


39~ 

10 ^ 17 
23 39 


62~ 

27 17 
62 39 


Table 16: ao 


= [2,3,2 


= 1/(2- 


f 1/(3 + 


1/(2 + r- 1 ))) ~ 0.4356 (360a ^ 


156.8°). 


a 


2 


3 

(2^5) 


5 
12 

(5,7) 


7 

17 

(5,12) 


9 

22 

(5,17) 


16 
39 

(17,22) 


25 
61 

(22,39) 


41 

100 

(39,61) 




(2^3) 


n c 
a 


3~ 
i i 

3 2 


5~ 

2 1 

5 2 


7~ 

2 3 
5 7 


12~ 

2^5 
5 12 


17~ 

2 7 
5 17 


22~ 

9^7 
22 17 


39~ 

9 16 
22 39 


61~ 

25 16 
61 39 



Table 17: a = [2,2,3] = 1/(2 + 1/(2 + 1/(3 + t" 1 ))) ~ 0.4100 (360a ^ 147.6°). 



Having prepared the mathematical relationship between the initial diver- 
gence ao, the final divergence a and the trace length n c , we are in a position 
to give an account of what is special about the golden angle. As shown 
below in Figs. [16] and [191 discontinuous change in phyllotactic fraction a, 
or phyllotactic transition, involves reconstruction of the vascular structure. 
Therefore, it is advantageous for a plant to suppress the transitions as few 
as possible. As internodes vary in length during growth, the trace length per 
internode n c may change accordingly. For instance, n c may depend on the 
plastochron ratio a ( Appendix A ). Patterns with a fraction that appears 
in many places of Fig. [7] are stable against occasional changes in n c . The 
lowest fraction that appears more than once is |. Thus, systems with initial 
divergence angle giving rise to stable fractions are most likely to survive. 
Among all possible values of ao, the initial divergence angle which suffers the 
least number of phyllotactic transitions is the golden angle ao = [2] (137.5 
degree s). Thi s is a summary of the evolutionary mechanism for the golden 
angle ( lOkabd ( 120111 )). Fig. [8] shows phyllotactic fractions resulting from var- 
ious representative values of ao while n c increases up to eleven. The number 
of phyllotactic transitions is indicated by a dashed line. In this example, 
initial divergence angles from 135° to 154° are most likely to be naturally 
selected. 
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Figure 8: Phyllotactic fractions resulting while n c increases to eleven are arranged verti- 
cally for eight representative values of the limit divergence ao . A dashed line is the number 
of phyllotactic transitions counted from a \ phyllotaxis. Initial divergences within a nar- 
row range around the golden angle ao = [2] ~ 0.382 (137.5°) are most likely to survive. 
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a 


2 


3 


5 


7 


12 


19 31 


50 


(?%) 


8 

(3,5) 


13 
(5,8) 


18 

(5,13) 


31 

(13,18) 


49 80 
(lo,ol) (ol,4yj 


129 
//in srA 


n c 


3 ~ 


Ft ~ 


8 rsji 


1 3 ~ 


1 8 ~ 


31 ~ 4Q ~ 


80 ~ 


a 


i i 


1 2 


3 2 


5 2 


5 7 


12 7 12 19 


31 19 


— r\j — 

3 2 


— r\j — 

3 5 


— r\j — 

8 5 


r^j — 

13 5 


13 ~ 18 


31 ~ 18 31 ~ 49 


80 ~ 49 


Table 18: a = 


[2,1,1,2] 


= 1/(2H 


-1/(1 + 1/(1 


+ 1/(2 + 


t- 1 )))) ~ 0.3876 (360a 


~ 139.5°). 


a 


1 


1 


1 


1 


2 


3 5 




(its) 


(M) 


6 

(1,5) 


(lie) 


13 

(6,7) 


19 32 

(6,13) (13,19) 


(19,32) 


n c 


3~ 


4~ 


5~ 


6~ 


7~ 


13~ 19~ 


32~ 




~ i 


0~i 


o~i 


~ i 

r 


2 1 

rsj — 

13 fi 


2^3 5^3 
13 ~ 19 32 ~ 19 


5 8 
— rNj — 

32 51 




Table 19: a = 


[6,2] = 


1/(6 + 1/(2 


+ r- 1 ) = 


0.1567 (360a ~ 56.4°) 





3. Precision of initial divergence angle 



As the evolutionary mechanism relies on statistical screening processes, 
it does not predict a limit divergence angle with unlimited precision. It is 
an empirical fact that divergence angles at the level of the shoot apex are 
regulated toward a mean value comparable with an i deal angle given by the 
formu la ([2]) after some trans ient fluctuations (IDaviesI (119391 ) : ISnow and Snow 
(119621 ); iBarabe et al.l (120101 )). Excepting initial fluctuations, the precision 
with which leaves are organized on the apical meristems is remarkable. It 
is undoubtedly controlled by genetics, though it may be slightly affected 



by lig ht stimuli depending on the orientation (IKumazawa and Kumazawa 



(11971I )). Twenty samples of young shoots of Erigeron sumatrensis (Suma- 
tran fleabane) show mean divergence angles from 137.23° to 137.97° wi th 
the sample average of 137.499 ± 0.212° ( Kumazawa and Kumazawal ( 1971 )). 
The mean divergence angle of the indi yidual plant may deviate statistical ly 
significantly from the ideal limit angle ( iMaksymowych and Ericksonl ( Il977l )). 
Sometimes there occur other ideal divergence angles than the normal golden 
angle of 137.5°. Phyllotaxis of Musa sa pientum ( banan a) changes with the 
age of the plant from | through | to | dSkutchl (Il927h ). This is consistent 
with a unique initial divergence of a = [2, 3] (Table El), which se e ms to 
be true for all species of Musa propagated vegetatively. iRutishauser (|l998h 
has presented a remarkably exotic pattern of Picea abies (Norway spruce) 
showing a (6, 13) phyllotaxis (Table 119]) . 
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The evolutionary mechanism predicts a correlation between the range of 
values of the initial divergence a® and the highest-order fraction a attained 
in evolutionary or phylogenetic processes. The correlation may seem strange 
at first glance, as it appears as an advanced correlation in developmental 
or ontogenetic processes of a plant; the precision of divergence angles on 
a young shoot is determined by the phyllotactic form at its maturity. This 
phylogenetic correlation is contrasted with the instantaneous correlation that 
causal models predict between the divergence angle and parastichies of the 
standing pattern. In general, divergence angles of a (3, 5) phyllotaxis are 
widely variable within | < a < §, whereas the range is narrowed to | < 
a < | when the parastichy pair is raised to (5,8). Remark that these are 
general results drawn from regularity of phyllotactic patterns. The ranges 
may be restricted further depending on specific assumptions of models. For 
instance, consider a regular pattern with a parastichy pair (1,2), which is 
realized for any value of divergence angle. According to Schwendener's model, 
however, (1,2) patterns for < «o < 0.36, i.e., from to 128.6 degrees, are 
not rea l ized, fo r a tran si tion t o a (2, 3) phyllotaxis i n terven es at «o = 0.36 
( jAdlerl ( jl974f ): iLevitovl ( jl99lh ; iDouadv and Couderl (jl996h ). See the top 
branch of the zigzag path in Fig.[2j The threshold angle «o = 0.36 specifically 
depends on geometrical assumptions, e.g. the circular shape of 'leaves' on the 
stem cylinder surface. Accordingly, the divergence angles for the parastichy 
pair (1,2) is predicted to vary continuously within 0.36 < a < ~, i.e., from 
128.6° to 180°. The range is narrowed substantially but still so wide that it is 
incompatible with observations that divergence angles are very close to 137.5° 
even in systems of low phyllotaxis. Causal models attain a target pattern 
with 137.5° by way of an almost opposite (1, 2) pattern with divergence of 
about 180°. Therefore, they cannot but allow the wide latitude of divergence 
angles for the (1, 2) pattern, in disagreement with precis e control of actua l 
systems (cf. Fig. 12(a) ). This is a very old problem which Ivan Itersonl (119071 ) 
(p. 247) was well aware of. Nonetheless, it has been left unnoticed despite 
a marginal rise of vario us causa l mod els in recent years. With reference to 
experimental evidence, I Church! (119041 ) (p. 340) remarks that already at a 
(2, 3) system the ideal angle is attained within an error of about one degree. 
The present model explains the non-correspondence between divergence angle 
cto and parastichy numbers by relating the allowed range of cto not with the 
parastichy numbers but with the highest order fraction a that the plant 
would attain in its mature state. 

By measuring initial divergence angles for thirty species of plants, iFujita 



29 



(119391 ) found that frequency distributions of the divergence angles are almost 
independent of the parastichy numbers. The divergence angles cluster in a 
narrow range. The width of the range quantifies the remarkable constancy 
of the diver gence angle ( Fig. 12(a) ). This results look puzzling from a causal 
viewpoint (IJeanl (119861 . I1994 VK By contrast, they are consistent with the 



evolutionary mechanism in that t he init i al an gle a® is independent of the 
parastichy numbers. According to Fujita ( 1939[). initial divergence angles for 
the main sequence fall within 138 ± 7° (IFuiital (Il939h ). irrespective of the 
parastichy pair. This corresponds to | < «o < § (135 to 144 degrees), which 
is as expected for the highest-order fraction of a = •jg (Table p. Similarly , an 
estimate of 99 ±4° for Cunninghamia lanceolata (China fir) (IFujital (119391 )) is 



consistent with \ < ao < I for a 



A in Table [2J and a narrow scattering of 



4 ^ ^ u ^ 7 _ n 
151 ±3° for (2, 5) phyllot axis at the apex of Cephalotaxus drupacea (Japanese 

plum yew) (IFuiital (Il937f ) ) is consistent with ^ < a < f for a = ^ (Table [6]). 

Let us make a general remark that parastichy does not substitute for di- 
vergence angle. The former depends on size and shape of the pattern unit or 
on a radial or internodal length scale. Therefore, several different parastichy 
pairs may be arbitrarily related to a single divergence angle. Parastichy num- 
bers given in Tables [TMT81 are the simplest pairs, which normally represent 
contact parastichies. 

Large fluctuations in the initial divergence ao may cause the phyllotac- 
tic transition in the vascular structure, even if the trace length n c is fixed 
constant. To suppress the transition that could happen, the divergence ao 
has to be restricted within one of the ranges determined by n c . For a fixed 
length of n c = 5, the fraction a is plotted against the initial divergence ao 
in Fig. |9j To maintain a | phyllotaxis, the initial divergence ao must stay 
| (from 120° to 144°); otherwise one would observe occa- 



within | 



< a < 



sional excursions to | (for ao < |) or | (for | < ao) in the midst of a steady 
course of the | phyllotaxis. Similarly, to maintain a phyllotaxis, the initial 



13 

divergence a has to be kept within | < a < § (from 135° to 144°); oth- 
erwise one would find ^ (for a < |) or y| (for | < a ) within the mature 
state of the ^ phyllotaxis (cf. Fig. |H]). Thus, it is explained why the ini- 
tial divergence angle has to be 'quantized' or fixed around a special constant 
with precision determined by the length of leaf traces. For this mechanism to 
work, stepwise changes in the fraction a of the vascular order, which are pre- 
sumed to occur if the initial divergence angle ao were not optimum, should 
incur penalties of extra energy. Thus, efficiency of the mechanism depends 
on the energy cost per transition, which should depend on species. By and 
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Figure 9: Phyllotactic fraction a versus initial divergence angle ag for a fixed length 
n c = 5 of leaf traces. The fractional order changes discontinuously while «o changes 
continuously. There are plateaus for five phyllotactic orders with |, |, |, | and |. The 
initial divergence «o is 'quantized' within a plateau to avoid the discontinuous transition. 
In other words, leaves are initiated regularly with a given angular precision. The wide 
plateau for a = | is the most unstable against changes in n c , while the plateau at a = | 
is the most stable. 



large, however, the number of transition may be used as a good measure of 
the total cost, at least first approximation (Fig. [8]). 



4. Fossil record and diversity of phyllotaxis 



Dickson! (1l87ll ) found that nine among thirteen specimens of fossil remains 



of Lepidodendron (scale tree) show helical phyllotaxis, of which only three be- 
long to the main sequence. This is in strikin g contrast to t he current domi- 
nance o f the main se quence in existing species (jFujital ( 119381 ) ; IZagorska-Marek 
( 1l985h : I.Teanl (jl994T )). Therefore, Dickson concluded that the phyllotaxis of 
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No. 



9,10 



11 



12 



13 



a 


13 

34 

(13,21) 


a 

55 

(21,34) 


bb 

144 

(55,89) 


21 

76 

(29,47) 


13 

60 

(23,37) 


21 
50 

(19,31) 


18 

59 

(23,36) 


47 

154 

(59,95) 




n c 
ato 


[21, 34) 

8^5 
21 13 

[2] 


[34,55) 

8 13 
21 34 

[2] 


[89,144) 

21 34 
55 89 

[2] 


[47,76) 

8 13 
29 47 

[3] 


[37,60) 

8^5 
37 23 

[4] 


[31,50) 

13 8 
31 19 

[2,2] 


[36,59) 

7 11 

23 36 

[3,3] 


[95,154) 

18 29 
59 95 

[3,3] 



Table 20: Ranges of n c and ao for the phyllotactic fract ion a and the c ontact parastichy 
pair (n,m) of the nine specimens of Lepidodendron bv iDicksonl (1871). Abbreviations 
[21,34) and |j- — ^ mean 21 < n c < 34 and ^- < a < jg. The bracket notation in @ 
is used for the limit divergence in the last row. Only the first three specimens belong to 
the main sequence ao = [2] (137.5°). 



Lepidodendron is extremely variable, as much so as that of those most vari- 
able plants like cacti. His results provide us with important information 
when they are analyzed in terms of the model. 

In the second and third line of Table [2U the phyllotactic fractions and 
the parastichy pairs for the nine specimens are presented after Dickson. The 
fourth and fifth line are the corresponding ranges of n c and ao according to 
Tables HJ |2j [3j |6l and |9j The last line is the limit divergence in terms of 
the bracket notation defined by (j2j) in the last section. For instance, in the 
second column, the specimen No. 1 has a || phyllotaxis (a = ||) with the 
parastichy pair of (13, 21), for which 21 < n c < 34 and ^ < ao < t|- The 
limit divergence of ao = [2] (137.5°) satisfies the latter condition. The spec- 
imens Nos. 1-3 belong to the main sequence ao — [2]- Fig. 10(a) represents 
graphically the parameter regions allowed for n c and ao- By comparison, 
Fig. 10(b) gives a theoretical result for the mos t favored regi ons in which the 
number of phyllotactic transitions is minimal (Okabe ( 201l[ )). 



According to Fig. 10(a) , the trace length n c appears to be independent of 
the initial divergence an. Moreover, n c is not as variable as ao- As the order 
of phyllotaxis is very high, there is considerable uncertainty in n c , while an 
is quite accurate. The specimens may be divided into two groups in terms 
of n c , i.e., one with n c ~ 50 and the other with n c > 100. The fact that the 
fossil specimens show various but accurate values of ao strongly suggests the 
evolutionary origin of the special divergence angles. It is impossible to tabu- 
late all phyllotactic fractions for such a large value as n c = 50 due to lack of 
space, but it is mentioned only that the number of possible phyllotactic frac- 
tions at n c = 50 amounts to 387 (~ "Sn^/it 2 /2). Among them, only the single 
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Figure 10: (a) Ranges of ao and n c for the specimens of lDicksonl (|l87lh (Table [20]) are 
painted black in the ao-n c plane, (b) A theoretical result for the regio ns in which t he 
number of phyllotactic transitions is minimal (adapted from Fig. 13 of lOkabel ( 20 1 ll ) ) ■ 
The golden angle ao — 0.382 (137.5°) is singled out for n c below Fibonacci numbers such 
as 34, 55, 89 and 144. 



fraction a = jk of the main sequence falls in the optimum regions depicted 
in Fig. 10(b) , while the specimens Nos. 8-12 do not meet the optimum condi- 
tion. Nevertheless, all the reported specimens possess the irrational numbers 
expressed in the form o f (j2J), as expected in the evolutionary mechanism (cf. 
Table 2 of lOkabel (|201lh ). Phyllotactic patterns for a = [3, 2], [2, 3], [2, 1, 2] 
and others are not reported, presumably because of lack of enough samples. 
Thus, anomalous patterns are regarded as relics of evolutionary processes. 

It has been an unresolved problem in what quantitative terms normal 
and anomalous phyllotaxis are differentiated. The number of phyllotactic 
transition during a steady growth provides us with a numerical measure of 
relative fitness in evolution. The most fit divergence angles a re ind icated in 
Fig. 10(b) They are pe aks of a 'fitness landscape' ( Niklas ( 1997f )). shown 
in Fig. [11] (Oka bd (201ll )). A close inspection of the frequency distribution 
curves o f iFuiital ( 11939m " indicates that a primary peak accompanies small sub- 
sidiary peaks at anomalous angles. In Fig. [TJ] Fujita's result for Lysimachia 
clethroides (gooseneck looses trife) is arrang ed along with transections of the 
fitness landscape in Fig. [11] iRobertsi ( 11984] ) has discussed that his chemical 
contact pressure model explains the anomalous subsidiary peaks. However, 
his conclusion is based on circular reasoning that anomalous systems are 
less frequent because they are anomalous. Similar fi tness c u rves are ob - 
tained for light absorption efficiency of rosette plants (INiklasI (119881 . Il998l ); 
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Figure 11: Three-dimensional fitness landscape. The vertical axis representing 'fitness' 
is n c minus the number of phyllotactic transition. The variables on the base plane are 
ao and n c . The number of transition increases with n c . The fitness has a flat bottom 
minimum in the worst case of ao — 0, whereas there are 'fitness peaks' at ao — [2] ~ 0.38 
(the main sequence), ao = [3] ~ 0.28 (an accessory sequence) and others, whose widths 
decrease as n c increases (|Okabd (|201lh i 
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(a) 



(b) 



Figure 12: (a) Frequency dist ributi o n of i nitial divergence angles for a (1,2) phyllotaxis 
of Lysimachia clethroides by Fuiita ( 19391 ). (b) An enlarged view of a relative 'fitness' 
in Fig. QT] is plotted against the divergence angle in degrees, 360a . The peak plateau 
extends from 135 to 144 degrees at n c — 12, and from 135 to 138 degrees at n c = 19. 



Pearcv and Yangl < jl99«h ; iKing et all ( 120041 ) ). 

Let us remark incredible precision of the divergence angle. As already 
mentioned, it is no less astonishing than the widely noticed fact that di- 
vergence angles converge on one of the special irrational numbers. Let us 
take the specimen No. 13 as an example. The divergence angle of the ^ 
phyllotaxis is a rational number 360a ~ 109.870°. This is very close to an 
irrational, ideal angle of a = [3,3], or 360a — 109.877°. According to 
Table [2UJ the range of a for the phyllotaxis is very narrow, that is, 



109.831° < 360a < 109.895°, 



(3) 



or 360a — 109.863 ± 0.032 degrees. The relative precision is less than about 
a part per three thousand. For reference, we present results that would be 
obtained if a happens to be off the narrow range of ([3]). Instead of ^ and 
the parastichy pair (59,95) for ©, we would have obtained ^ and (59,82) 
if cto were slightly below the lower limit of (j3J), or a — and (36, 95) if an 
were above the upper limit of ([3]). Neither of the last two cases is listed in 
Tables [TKITBl for they are hardly ever likely to occur. The plants' ability to 
distinguish ^ from ^ and ^ is due to high precision regulation of initial 
divergence angle. The range width of «o depends not so mu ch on an as on 
n c . Indeed, we find Aa ~ (r/n c ) 2 according to Eq. (B.39) in lOkabd ( 1201 it ). 
The precision as hi gh as the above ca n not be attained by a limited number 
of cells on the apex ( Koch et al.l ( 1998 ); Meinhardt et al.l ( 1998 ): Smith et al 
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(l2006a[ )). It seems very unlikely that existing causal models can explain this 
anomalous phyllotaxis with this precision in this probability of one out of 
thirteen specimens. 

Diversities of phyllotaxis is considered as a result of selective pressures be- 
ing ineffective. In extant plants, the main Fibonacci phyllot axis is dominant 
while some species specifically show very diverse phyllotaxis (jZagorska-Marek 
(119941 ) ) . In general, the trait diversity will be reduced if there is selective pres- 
sure acting on it. Strength of selective pressure depends specifically on extra 
cost required while rearranging phyllotactic patterns of leaf traces during 
growth of individual plants. Accordingly, the diversity may be preserved for 
some reason or other, e.g., when leaf traces are so fragile that the energy cost 
of rearrangement is insignificant. This view is consistent with recent research 
on Licopodium revealing a link b etween variability of leaf traces and diversity 
of phyllotaxis (jGola et al.l (120071 )). In contrast, the diversity in phyllotaxis of 
scale trees is considered as a result of strong selective pressure of insufficient 
time durations, strong because divergence angles are highly accurate whereas 
insufficient because various angles besides 137.5° are still in existence. In dis- 
cussing diversity of phyllotaxis, one should make a clear distinction between 
the variance, or standard deviation, of divergence angle of an individual and 
varieties of divergence angles of individuals. This section was devoted to the 
latter, while the former was discussed in the last section. 



5. Phyllotaxis and vascular organization 



Girolamil (119531 ) investigated the relation between phyllotaxis and vascu- 
lar organization of Linum (flax), whose vascular structures of a y|, ^ and 
phyllotaxis are given in the left, center and right of Fig. [131 respectively. 
On the one hand, the genetic spirals of the ^ and ^ phyllotaxis wind up to 
the right (counterclockwise), while it goes to the left (clockwise) for the ^ 
phyllotaxis. On the other hand, the main parastichies of the three patterns 
run in the same direction. That is to say, 5-parastichies for ^ (1-6-11-16-21, 
etc.), 7-parastichies for ^ (1-8-15-22-29, etc.) and 8-parastichies for ^ (1- 
8-15-22-29, etc.) run steeply from the bottom right to top left (clockwise). 
The most direct vascular connection goes along the main parastichies. The 
vascular bundles of these parastichies are recognized as sectioned clusters 
in a transverse section of the stem, called parastichy sectors. As shown in 
the inset of Fig. [131 the ^ phyllotaxis stem is divided int o five par a stichy 



sectors. In what follows, the following points remarked by I Girolamil (119531 ) 
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Figure 13: Diagrams of the primary vascular system of ^ (left), jg (center), ^- (right) 
phyllotaxis of Linum usitatissimum. In the inset (top left), dashed lines mark off f ive paras- 
tichy sectors in a transverse section of the phyllotaxis stem. Adapted from iGirolamil 
(ll95Sh . 
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are analyzed in terms of the model, whereby some useful general rules are 
pointed out: 

(Gl) The number of the parastichy sectors (5, 7 and 8 for y|, ^ and 
respectively) agrees with the numerator of the phyllotactic fraction for ^ 
and ^ of the main sequence, but not for ^ of the accessory sequence. 

(G2) The length of leaf traces per internode increases with the number of 
parastichy sectors, namely 12 for ^, 17 for ^, and 19 for Jy approximately 

(G3) There is no correlation in the relative directions of the genetic spiral 
and the parastichies. 

As noted in the first point, there is no easy-to-use general formula be- 
tween the parastichy numbers and the phyllotactic fraction (see below how- 
ever), but the numerical correspondence is immediately read from Fig. [7J 
and Tables [TMTHl According to Table [2, the parastichy pair of the ^ phyl- 
lotaxis is (7, 11). The number of parastichy sectors is the small number of the 
parastichy pair. Therefore, the number 7 of the parastichy sectors of the j| 
phyllotaxis is obtained. Unlike the numerator, the denominator satisfies sim- 
ple rules. Most notably, the denominator of a fraction is equal to the sum of 
the contact parastichy pair corresponding to the fraction (e.g. 18 = 7 + 11). 
Mathematical relations between various numbers in phyllotaxis have been 
investigated since early times on an e mpirical ground based on pu r ely mathe- 



matic al prop e rties o f a regular lattice ( iBravais and Bravaisl (118371 1 ; iNaumann 
f |l845l ): l,Teanl (|l994h ). 

On the second point, Tables CD and [2] give the conditions 8 < n c < 13, 
11 < n c < 18 and 13 < n c < 21 for the phyllotactic fractions y| and 
Jy, respectively. The predictions of the model are supported by the reported 
values n c = 12, 17 and 19, which satisfy their respective conditions near 
their upper limits. Nevertheless, a close look at Fig. [T3] indicates that these 
figures are not accurate. As a matter of fact, n c appears not constant but 
somewhat larger in the upper part of the stem. Chang es in length of t he lea f 
traces are revealed in a more sophisticated analysis of iMeicenheimerl (119861 ). 
where progressive transitions from | through | and | up to ^ have been 
reported. Phyllotactic transition caused by changes in n c is discussed in the 
next section. 

On the third point, a general rule holding between directions of paras- 
tichies and the genetic spiral is presented based on Fig. [71 To this end, 
it is convenient to introduce a 'mother' fraction of a fraction a, which is 
defined as the fraction lying immediately below the fraction a in the tree 
of Fig. [7J The mother fractions of A and ^ are |, ^ and ^, respec- 
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Figure 14: Three branches from Fig.[7]with which to explain spiral directions of j*, ^g, ^ 
phyllotaxis. As the fraction ^ and are numerically bigger than their 'mother' fraction 
YY and | , their main parastichies of 7 and 5 are contrary in direction to the genetic spiral. 
On the contrary, 8 parastichies for ^j- are in the same direction as the genetic spiral. 



tively. It is shown that if and only if a phyllotactic fraction a is numerically 
bigger than its mother fraction, the main parastichies run in the direction 
opposite to the genetic spiral. (The main parastichies are gentle, long spi- 
rals characterized by the small number of the contact parastichy pair.) The 
fraction a = y| and ^ are bigger than the mother fraction | and A, respec- 
tively, while a = ^ is smaller than the mother fraction ^. The magnitude 
relations are schematically shown in Fig. |5] extracted from Fig. [7J Thus, 
the above rule explains Girolami's observation consistently. In practice, this 
rule may be used to identify the direction of the genetic spiral of a high 
order phyllotactic pattern for which parastichies are far easy to fo llow vi- 
sually . Some s pecial cases of this general ru le have been remarked (I Church 
( 119041 ) (p. 9 6). iNamboodiri and Beckl (119681)') and occa sionally taken up for 
discussion (IMeicenheimerl ( 119861 ); iFredeen et al.l ( 120021 )). The directional re- 
lations between various spirals of a phyllotactic pattern are also mathematical 
consequences of the regularity of the phyllotactic pattern. 

The mother fraction enables us to state general rules for the phyllotactic 
fraction and the parastichy number: One of the parastichy pair for a frac- 
tion a is equal to the denominator of the mother fraction of a; The other 
number in the pair is determined such that the sum of the pair is equal to 



the denominator of a. Consider a 



18' 



for instance. One of its parastichy 
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Figure 15: When the final divergence a is numerically bigger than the initial divergence 
ao, the stem is twisted in the direction of the genetic spiral. N.B. Divergence angles are 
less than 180 degrees. 



pair is the denominator 11 of the mother fraction while the other is the 
difference of the denominators, 18 — 11 = 7. As a result, the parastichy pair 
(7, 11) is obtained for ^. Thus, the rules are used to relate the parastichy 
numbers and the phyllotactic fraction. 

The vascular systems shown in Fig. [TS] form closed networks. In each 
system, connections between leaf traces are formed along both the paired 
parastichies, so that the vascular bundles are divided into parastichy sec- 
tors. Among dicotyledons wi th helical phyllot axis, however, an open vascu- 
lar system is rather common (IBeck et al.l ( 119821 ) ) . Primit i ve ang iosperms and 
m any Kymnosperms have open vascular systems (IBeck! (120101 ) ) . According 
to IBeck et al.l ( 119821 ). open systems of five sympodia (a | phyllotaxis) char- 



acterize 67% of the species with helical phyllotaxy and are clearly a common 
type among dicotyledons. In an open system, leaf traces are connected along 
one direction. Although the present model determines the basic architecture 
of vascular phyllotaxis, it does not specify detailed structure of the reticulate 
pattern, whether it remains open or becomes closed. This is not a shortcom- 
ing of the model, because actual linkages betw een leaf traces are likely to be 
secondary events depending on circumstances (IKang et al. rfiooi n. 

To conclude this section, let us remark another obvious correlation be- 
tween the direction of the genetic spiral and the secondary torsion of the 
stem. The initial divergence «o is related to the fractional divergence a of a 
mature pattern by the angle of twist a — ao undergone in the secondary tor- 
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sion. The direction of the torsion is the same as the genetic spiral if and only 
if a > otQ. This is shown schematically in Fig. The direction of the sec- 
ondary torsion would not be difficult to check experimentally. In most typical 
cases, the direction is reversed, or the sign of a — a n chan g es, as n c crosses 
a threshold of phyllotactic transition. iBravais and Bravaisl (118371 ) evaluated 
the limit divergence ao from mature shoots by correcting the torsion angle 
a — otQ. 



6. Phyllotactic transition 

Larson ( 19771 ) has investigated phyllotactic transition in the vascular sys- 



tem of Populus (cottonwood). His result showing transition from a | to 
phyllotaxis is reproduced in Fig. (THl Each leaf has three traces; central, right 
and left traces are indicated with crosses, filled and open triangles, respec- 
tively. The leaf traces are connected with the stem vascular bundles to make 
sympodia. The sympodia are separated from each other, or the vascular 
system is open. The three traces leading to each leaf primordium arise on 
different sympodia. The number of the sympodia changes from five in the 
lower portion to eight in the upper portion of Fig. [16j The number agrees 
with the denominator of the phyllotactic fraction in each part. The region of 
the | phyllotaxis occurs in the basal stem above some primary leaves, while 
the | phyllotaxis occurs at mid and up per stem levels, principally in the 
zone of expanding leaves (jLarson (Il977h . In Fig.dSl once the transition is 
initiated at a point IA on a sympodium number 2, it progresses through the 
sympodia at points IB through IE. Three new central traces to establish the 
three additional sympodia of the | system are derived from left traces in se- 
quence at points IIA-IIC. Variou s interre l ation s between phyllotaxis and leaf 
development have been studied (jLarsonl (119801 ) ) . In what follows, a correla- 
tion between phyllotactic transition and lengths of the leaf traces is analyzed 
by means of the model, whereby supporting evidence of the model is pointed 
out. 

The lengths per internode of the leaf traces are optically read from Fig. [16] 
and plotted in Fig. [171 Arrows indicating the transition region between 
the | and I phyllo taxis in Fig. [T7] are marked in accordance with Fig. [TBI 



after iLarsonl (119771 ). By comparison, a dashed line at n c = 5 is drawn to 
indicate the theoretical threshold between the | and | phyllotaxis (Table 
[Q). In accordance with the model, the phyllotactic transition is triggered 
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SYMPODIUM NUMBER 3/8 PHYLLOTAXY 




SYMPODIUM NUMBER 2/5 PHYLLOTAXY 



Figure 16: Transition in the primary vascular system of a cottonwood plant from a | 
to | phyllotaxis. Central, ri ght and left traces are indicated with crosses, filled and open 
triangles, respectively. After iLarsonl (|l977t ). 
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Figure 17: Length per internode of the leaf traces in Fig. [16] is plotted against the leaf 
node number, the vertical axis of Fig. 1161 Arrows for a transition region between the | 
and | phyllotaxis and labels IA-IE and IIA-IIC to indicate initiation of the transition are 
marked in accordance with Fig. [16] by Larson. The phyllotactic transition is consistent 
with the threshold value of n c ~ 5 predicted by the model (Table [lj . 



by the increasing length of the leaf traces crossing a threshold value of five 
internodes. 

According to Table CD, phyllotactic transition is predictable. Transitions 
of the main sequence occur whenever the trace length n c crosses Fibonacci 
numbers. The trace length, li ke other param eters of the plant, is predictably 
correlated with plant vigor (ILarsonl ( 1l980l )). Therefore, in principle, the 
model allows us to control phyllotaxis artificially. In Sec. [2j leaf traces are 
assumed to have a common length. As noted at the end of the last section, 
the direction of the secondary torsion is reversed when n c crosses a threshold 
value, so that it may be fixed by a leaf trace of length longer than the 
threshold. Fig. [T8l schematically shows that long leaf traces 10, 11 and 12 
trigger a transition from | to |. In the transition region of Fig. [T71 three left 
traces (open triangles) of the node number 7, 8 and 9 are the first to cross the 
threshold at n c = 5. These are the very traces labeled with IIA, IIB and IIC 
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Figure 18: A phyllotactic pattern with aa = 1/(1 + r) (cf. Fig. 3(a)| . Length of leaf 
traces is n c = 4 (solid bars) except for 10, 11 and 12 with n c = 6 (bold bars). The longer 
traces can induce a transition from a = | in the lower portion (cf. Fig. |3(b)| to a = | 
in the upper portion (cf. Fig. 5(b)). At the transition, the longer traces deflect main 



parastichies (dotted lines), and the parastichy number increases from 5 to 8. 



by Larson as those from which the three extra sympodia branch. A close look 
at Fig. [in] reveals that central traces below and above the transition region 
are inclined in the opposite direction. This is consistent with the prediction 
of the model, for a — | < «o < §■ Furthermore, five right traces (filled 
triangles) striking around n c ~ 6 in Fig. [T7] agree with the special traces 
labeled with IA through IE. Thus, the observation supports the special role 
of the Fibonacci number 5 for the trace length n c . 

Two-step transition from a | to & phyllotaxis is shown in Fig. IT9l after 



Larson! (119771 ). where steady increase in length of leaf traces is more obvious 
than Fig. [16] Fig. [20] is obtained from Fig. [19] in the same manner as Fig. [T7] 
is obtained from Fig. HE] Leaf positions at which the I phyllo taxis starts 
and ends are marked on the right side of Fig. IT91 by iLarsonl (119771 ). according 
to which the transient pattern of the | phyllotaxis is maintained for the 
leaves with plastochron index from 5 to —7. Accordingly, the corresponding 
positions are marked by arrows in Fig. [20] On the other hand, horizontal 
lines at Fibonacci numbers 3, 5, and 8 in Fig. [20] theoretically divide the 



regions for the |, |, | and ^ phyllotaxis (Tabled]). Thus, it is confirmed 
again that continuous changes in length of leaf traces cause discontinuous 
transitions in the vascular structure. 

Fig. [20] indicates that n c increases steadily up to an upper bound of about 



44 




45 





13 




12 




1 1 


LU 


1 


o 




RA 


9 


h- 




LL 


8 


< 




LU 




I 


7 


LL 




o 




X 


6 


h- 




ON 


5 


LU 






4 




3 




2 




1 



i r 



~i r 

3/8 start 



~i 1 r 

3/8 end 



5/13 PHYLLOTAXIS 

3/8 

*. .. A..A.... A A .. A 

A A A A 

2/5 a 

A 

-x-x-x-^x-x-txr^^-x-g- 



A AAA' 



A a A A A ■ 



XX a .XaX*xX X X" 



x x *x 

AAA* 



_1_ 



_1_ 



12 



4 -4 -8-12 

LEAF PLASTOCHRON INDEX 



-16 



-20 



-24 



Figure 20: Length per internode of the leaf traces read from Fig. [19] is plotted against 
the leaf index (the vertical axis of Fig. [T 9|. Two arrow s at the top indicate where the | 
phyllotaxis starts and ends according to Larson ( 19771 ) (see Fig. Qj|]). According to the 



theoretical model, stable regions for the | , | and phyllotaxis are separated by horizontal 
dashed lines at Fibonacci numbers 3, 5 and 8 (cf. Table [TJ. Thus, Larson's estimate of the 
region of the | phyllotaxis agrees with the theory. Left traces (open triangles) reaching a 
maximum length of about 10 internodes is consistent with an observation that the highest- 
order phyllotactic fraction that this plant attains is 



10. This observation is consistent with the fact that the ^ phyllotaxis was 
the stable pattern of the old plant ( Larson ( 1980l )). According to the model, 
the y| phyllotaxis is stable insofar as n c lies between 8 and 13, i.e., there 
is a 5-internode allowance for the trace length of the ^ phyllotaxis. The 
main sequence is special for this wide cleara nce between s uccessive threshold 
values. The interval is denoted as An c in lOkabd (1201 ll ). As shown there, 
the widest clearances are achieved for Fibonacci numbers, and a sequence of 
Fibonacci numbers is realized when the limit divergence angle is one of the 
special irrational numbers related to the golden ratio. As shown in Fig. [HJ 
the number of transitions encountered while n c grows up to above 10 is kept 
to a minimum number insofar as the initial divergence is restricted within 
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I < a < | (from 135° to 154°, as noted at the end of Sec. EJ). When 
n c becomes larger than 12, the range is narrowed to | < «o < § (from 
135° to 144°). Thus, the normal phyllotaxis of the main sequence is singled 
out. Owing to the observation that the highest-order fraction was a = j^, 
the model predicts that the initial divergence an should be contained within 
| < a < |, just as observed by iFujital (119391 ) for other species (Sec. |2J). 
Unfortunately, initial divergence s of the co t tonwq od plant are not available 
to us. To support this argument, iPulawskal ( 119651 ) has reported for Actinidia 
arguta (hardy kiwi) that initial divergence remains constant despite changes 
in the vascular organization between |, ^ and ^. 

When n c is increased past 8, the model predicts vascular phyllotaxis of 
either a = | or q = | depending on whether the initial divergence cto is 
smaller or larger than | (angle of 144°). Suppose «o = [3] (99.5°), then one 
should have five threshold lines at 2, 3, 4, 7 and 11 (Table |2J), instead of 
four thresholds at 2, 3, 5 and 8 for cto = [2] in Fig. [201 If the initial diver- 
gence were «o = [5] (64.1° in Table HJ), one should have six threshold lines 
at n c = 2, 3, 4, 5, 6 and 11 separating patterns of a = |, |, |, |, |, ^ and 
^= (cf. Fig. [8]). The vascular phyllotaxis is very unstable. The instability is 
energetically unfavorable. Therefore, «o = [5] (64.1°) is very improbable to 
survive natural selection because of the multiplicity of expected transitions. 
A general remark should be made when discussing multiple patterns in se- 
quence. In order for a pattern with a definite value of a to be distinguished 
as such, the pattern should consist of more leaves than the denominator of 
the fraction a. This holds true if n c varies sufficiently gradually; otherwise 
phyllotaxis transition may not be distinctly discernible. 

Last but not least, whorled phyllotaxis has not been discussed in this pa- 
per. A J-jugate pattern with J fundamental spirals is formed when J leaves 
are borne at each node. Compared with a helical phyllotaxis with J = 1, di- 
vergence angles of a J-jugate system are divided by J and the parastichy pairs 
(m, n) are multiplied by J. Therefore, one obtains < J«o < \ an d J(m, n) 
for the divergence angle and parastichy pair of a J-jugate system. It is known 
that sometimes vascular structure may change between helical and whorled 
phyllotaxis during ontogeny. This type of 'anomalous' phyllot actic transitio n 
also appears to be caused by a de c rease in length of leaf traces (jJensenl (119681 ); 
Beck et al.l (119821 ); iKwiatkowskal (119951 )). The present model gives a — \ for 
1 < n c < 2 and a = | for 2 < n c < 3 irrespective of a . Correspond- 
ingly, it seems natural to consider that a whorled phyllotaxis is a variation 
of the most primitive alternate phyllotaxis and that a whorled phyllotaxis is 
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triggered as n c becomes less than 1. However, changes in the vascular struc- 
ture have to be coordinated with changes in the po sitioning of initia t ed lea f 
primordia whi l e a whorled pattern is e stabli shed (jZagorska-Marekl ( 11994] ); 
Meicenheimen dl998h : iKellv and Cooke! fl2003h l The physiological processes 
involved are unlikely to be amenable to simple mathematical analysis. Still, 
a similar transition rule as a helical pattern should hold for an established 
whorled pattern in terms of trace length redefined with a new internode. 



7. Conclusions 

The present work puts forward an important role of Fibonacci numbers 
as critical values of the length per internode of leaf traces played in vascular 
phyllotaxis transition. 

The regular arrangement of leaves and the regularity in divergence angle 
of 137.5° are a result of selective pressure to reduce possible changes in the 
vascular structure during growth, i.e., aperiodic arrangements will necessitate 
extra nutrients to reconstruct the sectorial or fractional order of vascular 
connections. 

The phyllotactic fraction a of mature patterns of leaf traces normally 
makes transitions through |, |, |, |, ^, whenever the number of 

internodes traversed by the leaf traces, n c , crosses Fibonacci numbers, 1, 
2, 3, 5, 8, 13, 21, • • • . The Fibonacci numbers make appearances because 
initial divergence angle «o of leaves at the shoot apex is normally the golden 
angle of about 137.5° with a good precision. The golden angle is prevalent 
because it is the selectively advantageous angle at which the number of the 
phyllotactic transition is the minimum (Fig. [8]). The precision of the initial 
divergence is determined by the trace length n c . 
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Appendix A. Relation between the trace length n c and the plas- 
tochron ratio a 

A point on a cylinder surface is located with the angular coordinate <p 
and the height z. Leaves on a stem are represented by a lattice of points 
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given by ip = 2iran (in radians) and z = hn, where a is a constant angle 
of divergence, h an internode length, and n is an integer index. On the 
other hand, a point on a plane is located in a polar coordinate system (r, <p), 
where r and (p is the radial distance from the central axis and the angular 
coordinate about the axis, respectively Leaf primordia at a shoot apex are 
represented by r = a n and (p = 2Tcan, where a is a plastochron ratio. In a 
conformal growth preserving angles, the two representations are related by 
2txz = logr. Hence, the internode length h corresponds to the logarithm of 
the plastochron ratio, ^loga. The number of internodes traversed by the 
leaf traces is n c = Z\ t /h = 2nZ\ t /\oga 1 where Z\ t is a length of leaf traces in 
the stem. Therefore, n c may be regarded as inversely proportional to log a, 
or the relative growth rate per plastochron The growth rate should 

depend on cell types. Accordingly, n c may change during plant growth. 
The plastochron r atio may chang e as a result of alteration in size of the 



apex and primordia. iRichardsl ( I195ll ) discussed changing phyllotaxis to the 



effect that a continuous shift in the parastichy pair of normal Fibonacci phyl- 
lotaxis is linearly correlated with a double logarithm log(loga). He defined 
the phyllotaxis index (P.I.) by 

P.I. = 0.38 - 2.39 log 10 log 10 a, (A. 1) 

where numerical values are chosen such that the index assumes an integral 
value whenever two sets of parastichies in the Fibonacci system intersect 
orthogonally. The crossing angle between the contact parastichies changes 
continuously as a function of the plastochron ratio. In this descriptive model, 
the divergence angle ao is fixed at the golden angle. 

Changing phyllotaxis due to change in the plastochron ratio is consistent 
with the present model of vascular phyllotaxis. In this model, the divergence 
a on the stem changes discontinuously, however. To show a correspondence 
between changes in phyllotaxis on the apex and the stem, let us consider the 
normal phyllotaxis with an initial divergence of the golden angle ao = t~ 2 
(Tabled]). Let us introduce the Fibonacci sequence F n generated from initial 
integers F\ = 1 and F 2 = 1 by the recurrence relation F n+2 = F n+1 + 
F n . Accordingly, F n = 1,1,2,3,5,8 and 13 for n = 1,2,3,4,5,6 and 7, 
respectively. In terms of F n , the phyllotactic fraction a = and the 

parastichy pair (F n , F n+ i) are obtained for F n+ i < n c < F n+2 , or for 

{n + 1) log r — log \/5 < log n c < (n + 2) log r — log y/h 
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owing to an approxima te formula F n ~ r n /\/5 valid for large n (see below 



(B.31) in lOkabd ( l201lh ). Therefore, logn c is proportional to the integer 
index n. 

To put it concretely, we get a = | and the parastichy pair (1,2) for 
2 < n c < 3 or 

0.7 < logn c < 1.1, 
a = | and (2, 3) for 3 < n c < 5 or 

1.1 < logn c < 1.6, 
a = | and (3, 5) for 5 < n c < 8 or 

1.6 < logn c < 2.1, 
a = an d (5, 8) for 8 < n c < 13 or 

2.1 < logn c < 2.6, 

and so on. Thus, the shift in the parastichy pair is linearly correlated with 
logra c oc log(loga). This is a general property holding also for other initial 
divergences found in nature. 

For the systematic study of the mature stem, the index n c is more usefully 
regarded as a developmental index than a, not only because an internode is 
a natural unit of length as the plastochron is the developmental unit of time, 
but values of n c allowed for a phyllotactic pattern are delimited by the spe- 
cial integers traditionally familiar to those who are enchanted by phyllotaxis; 
Fibonacci numbers. For a given initial divergence, the numbers comprise a 
sequence generated by the Fibonacci recurrence relation F n+2 = F n+ \ + F n 
from a pair of different seed integers. The main sequence, 1, 2, 3, 5, 8, • • • 
in Table [U is generated from the simplest seed pair (1, 2). The next simplest 
seed integers (1,3) give the accessory sequence 1, 3, 4, 7, 11, 18 • ■ ■ of Ta- 
ble |2j In this manner, any phyllotactic sequence is characterized by a pair 
of seed integers, as well as the limit divergence o:q. This is in accordance 
with accumulated empirical wisdom of phyllotaxis. Traditionally, these spe- 
cial integers have been remarked in connection with parastichy numbers (cf. 
Tables IT HTB"]) . The present work puts emphasis on these numbers as critical 
values for the length per internode of leaf traces. This point has never been 
remarked before. 
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